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ABSTRACT 

The project Massive Unseen Companions to Hot Faint Underluminous Stars from SDSS (MUCHFUSS) aims at finding sdBs with 
compact companions like supermassive white dwarfs (M > 1.0 M G ), neutron stars or black holes. The existence of such systems is 
predicted by binary evolution theory and recent discoveries indicate that they are likely to exist in our Galaxy. 

A determination of the orbital parameters is sufficient to put a lower limit on the companion mass by calculating the binary mass 
function. If this lower limit exceeds the Chandrasekhar mass and no sign of a companion is visible in the spectra, the existence of 
a massive compact companion is proven without the need for any additional assumptions. We identified about 1100 hot subdwarf 
stars from the SDSS by colour selection and visual inspection of their spectra. Stars with high velocities have been reobserved and 
individual SDSS spectra have been analysed. In total 127 radial velocity variable subdwarfs have been discovered. Binaries with 
high RV shifts and binaries with moderate shifts within short timespans have the highest probability of hosting massive compact 
companions. Atmospheric parameters of 69 hot subdwarfs in these binary systems have been determined by means of a quantitative 
spectral analysis. The atmospheric parameter distribution of the selected sample does not differ from previously studied samples of 
hot subdwarfs. The systems are considered the best candidates to search for massive compact companions by follow-up time resolved 
spectroscopy. 
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! 1 ■ Introduction to ra nge from 39 % to 78 % (e.g. Maxted et al. l2001l Napiwotzki 

■ et al. 1200 4a). Several studies were undertaken to determine the 

■ Subuminous B stars (sdBs) are core helium-burning stars with orbital parame ters of subdwarf binaries (e.g. Edelmann et al. 
; veryjhin hydrogen envelopes and masses around 0.5 M (Heber 2005 . M orales-Rueda et al. 2003a). The orbital periods range 

L; , |19861l. A large fraction of the sdB stars are membe rs of sh ort pe- from o 07 to > 10 d with a peak at 5 - 1 d 
. £h ' nod binaries (Maxted et al. 120011 Napiwotzki et al. 2004ai. After 

^ I the discovery of close binary subdwarfs, several studies aimed at For c i ose binary sdBs common env elope ejectio n is the most 

H ' determining the fraction of hot subdwarfs residing in such sys- pro bable formation channel (Han et al. l2TjQ2l 125531 In this sce- 

. terns. Samples of hot subdwarfs have been checked for radial ve- nar i two ma i n sequence stars of different masses evolve in a 

locity (RV) variations. The binary fraction has been determined binary systerrl . The more massive one will reach the red giant 

phase first and fill its Roche lobe near the tip of the red-giant 

Send offprint requests to: S. Geier, branch. If the mass transfer to the companion is dynamically un- 
e-mail: geier@sternwarte.uni-erlangen.de stable, a common envelope is formed. Due to friction the two 
* Based on observations at the Paranal Observatory of the European ste Har cores lose orbital energy, which is deposited within the en- 
Southern Observatory for programme number 081.D-0819. Based on velope and leads to a shortening of the binary period. Eventually 
observations at the La Silla Observatory of the European Southern ^ common envelope is ejected and a close binary system is 
Observatory for programme number 082.D-0649. Based on observa- £ , , . . , ,. , . JT -, , 

' , e - . * . , . tt' ai - att ax . formed, which contains a core helium-burning sdB and a main 
tions collected at the Centra Astronomico Hispano Aleman (CAHA) at . , . . . ° 

Calar Alto, operated jointly by the Max-Planck Institut fur Astronomie sequence companion. A binary consisting of a main sequence 

and the Institute de Astroflsica de Andalucfa (CSIC). Based on ob- star and a whlte dwarf ma y evolve t0 a close blnar y sdB wlth a 

servations with the William Herschel Telescope operated by the Isaac white dwarf companion in a similar way. Only in very special 

Newton Group at the Observatorio del Roque de los Muchachos of the and hence rare cases tight constraints can be put on the nature 

Instituto de Astrofisica de Canarias on the island of La Palma, Spain. of the companions, that is if the systems are eclipsing or show 
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Fig. 1. The RV semiamplitudes of all known sdB binaries with 
spectroscopic solutions plotted against their orbital periods (see 
Table lA.lt . Binaries which have initially been discovered in pho- 
tometric surveys due to indicative features in their light curves 
(eclipses, reflection effects, ellipsoidal variations) are marked 
with open circles. Binaries discovered by detection of RV varia- 
tion from time resolved spectroscopy are marked with filled di- 
amonds. The dashed, dotted and solid lines mark the regions to 
the right where the minimum companion masses derived from 
the binary mass function (assuming 0.47 M Q for the sdBs) ex- 
ceed 0.45 Mq, l.OOM and 1.40M o . The two post-RGB objects 
in the sample have been excluded, because their primary masses 
are much lower. 



Due to the tidal influence of the companion in close binary 
systems, the rotation of the primarjQ becomes synchronised to 
its orbital motion. In this case it is possible to constrain the mass 
of the companion, if mass, projected rotational velocity and sur- 
face gravity of the sdB are known. Geier et al. (2008. 2010a. 
2010b) analysed high resolution spectra of 41 sdB stars in close 
binaries, half of all systems with known orbital parameters. In 
3 1 cases, the mass and nature of the unseen companions could 
be constrained. While most of the derived companion masses 
were consistent with either late main sequence stars or white 
dwarfs, the compact companions of some sdBs may be either 
massive white dwarfs, neutron stars (NS) or stellar mass black 
holes (BH). However, Geier et al. d2010bl l also showed that the 
assumption of orbital synchronisation in close sdB binaries is not 
always justified and that their sample suffers from huge selection 
effects. 

The existence of sdB+NS/BH systems is predicted by bi- 
nary evolution theory (Podsiadlowski et al. 120021 Pfahl et al. 
2003). The formation channel includes two phases of unstable 
mass transfer and one supernova explosion. The predicted frac- 
tion of sdB+NS/BH systems ranges from about 1% to 2% of the 
close sdB binaries (Geier et al. 1201 Obi Yungelson & Tutukov 
120531 Nelemans l20T0l i. 

2. Project overview 

The work of Geier et al. (12010b) indicates that a population of 
non-interacting binaries with massive compact companions may 
be present in our Galaxy. The candidate sdB+NS/BH binaries 
have low orbital inclinations (15 - 30°, Geier et al. 2010b ). High 
inclination systems must exist as well. A determination of the or- 
bital parameters allows one to put a lower limit to the companion 
mass by calculating the binary mass function. 



(M comp + M sdB ) 2 2nG 



other indicative features in their light curves (see the catalogue 
of Ritter & Kolb 2009 and references therein). 

Subdwarf binaries with massive WD companions turned out 
to be candidates for SN la progenitors because these systems 
lose angular momentum due to the emission of gravitational 
waves and start mass transfer. The mass transfer, or the subse- 
quent merger of the system, may cause the WD to approach the 
Chandrasekhar limit, ignite carbon under degenerate conditions 
and explode as a SN la (Webbink [19841 Iben & Tutukov [19841) . 
One of the best known candidate system for this double degen- 
erate merger scenario is the sdB+WD binary KPD 1930+2752 
(Maxted et al. l2000al Geier et al. l2007l ). Mereghetti et al. ([20091 
showed that in the X-ray binary HD 49798 a massive (> 1 .2 M ) 
white dwarf accretes matter from a closely orbiting subdwarf O 
companion. The predicted amount of accreted material is suffi- 
cient for the WD to reach the Chandrasekhar limit. This makes 
HD 49798 another candidate SN la progenitor, should the com- 
panion be a C/O white dwarf (Wang et al. 2009). SN la play a 
key role in the study of cosmic evolution. They are utilised as 
standard candles for determining the cosmological parameters 
(e.g. Riess et al. [19981 Leib undgu t llOOTl Perlmutter et al. [T999l) . 
Most recently Perets et al. (120101 1 showed that helium accretion 
onto a white dwarf may be responsible for a subclass of faint and 
calcium-rich SN lb events. 



The RV semi-amplitude K and the period P can be derived 
from the RV curve; the sdB mass M s dB, the companion mass 
M comp and the inclination angle i remain free parameters. We 
adopt M sd B = 0.47 M Q and i < 90° to derive a lower limit for 
the companion mass. Depending on this minimum mass a quali- 
tative classification of the companions' nature is possible in cer- 
tain cases. For minimum companion masses lower than 0.45 M 
a main sequence companion can not be excluded because its lu- 
minosity would be too low to be detectable in the spectra (Lisker 
et al. 120051 1. If the minimum companion mass exceeds 0.45 M 
and no spectral signatures of the companion are visible, it must 
be a compact object. If it exceeds the Chandrasekhar mass and 
no sign of a companion is visible in the spectra, the existence of a 
massive compact companion is proven without the need for any 
additional assumptions. This is possible, if such a binary is seen 
at high inclination. The project Massive Unseen Companions 
to Hot Faint Underluminous Stars from SDSS0 (MUCHFUSS) 
aims at finding sdBs with compact companions like supermas- 
sive white dwarfs (M > 1.0 M ), neutron stars or black holes. 



1 The more massive component of a binary is usually defined as the 
primary. However, in most close sdB binaries with unseen companions 
the masses are unknown and it is not possible to decide a priori which 
component is the most massive one. For this reason we call the visible 
sdB component of the binaries the primary throughout this paper. 

2 Sloan Digital Sky Survey 
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Fig. 2. Left panel. SDSS g - r-colours plotted against u - g of all stars. The grey dots mark all stellar objects with spectra available in 
the SDSS database. Most of them are classified as DA white dwarfs. The solid diamonds mark (He-)sdO stars, the solid squares sdB 
and sdOB stars. Open squares mark hot subdwarfs with main sequence companions visible in the spectra. Most of these objects are 
white dwarfs of DA type. Right panel. Only subdwarfs with g < 18 mag are plotted. The sequence of composite objects is clearly 
separated from the single-lined stars. Synthetic colours from Castelli & Kurucz (2003) for stars with temperatures ranging from 
14 000 K to 50 000 K (logg = 5.0) are marked with upward triangles and connected. The stepsize of the colour grid is 1000 K. The 
labels mark models of certain temperatures. 



First results of our follow-up campaign are published in Geier et 
al. d20TT1 ). 

There is an interesting spin-off from this project: The same 
selection criteria that we applied to find such binaries are also 
well suited to single out hot subdwarf stars with constant high 
radial velocities in the Galactic halo like extreme population II 
stars or even hypervelocity stars. To refer to this aspect we coin 
the term Hyper-MUCHFUSS for the extended project. First re- 
sults are presented in Tillich et al. (1201 1) . 



3. Target selection 

The high fraction of sdB stars in close binary systems was ini- 
tially discovered by the detection of RV shifts using time re- 
solved spectroscopy (Maxted et al. 1200 11 1. In the past decade 
about 80 of these systems have been reobserved and their orbital 
parameters determined. We summarize the orbital parameters of 
all known sdB binaries and give references in Table lA.fl (see also 
Fig.0. 

As far as the companion masses of the known sdB binaries 
could be constrained, it turned out that most companions should 
be either late main sequence stars with masses lower than half 
a solar mass or compact objects like white dwarfs. Targets for 
spectroscopic follow-up were selected in different ways depen- 
dent on the specific aims of each project. 

For the MUCHFUSS project the target selection is optimised 
to find massive compact companions in close orbits around sdB 
stars. In order to discover rare objects applying the selection cri- 
teria explained in the forthcoming sections, a huge initial dataset 
is necessary. The enormous SDSS database (Data Release 6, 
DR6) is therefore the starting point for our survey. Best sky 



Table 1. Survey observations. The first column lists the date of 
observation, while in the second the used telescope and instru- 
mentation is shown. In the third column the initials of the ob- 
servers are given. 



Date 



Telescope & Instrument Observers 



January- June 2008 
2008/04/29-2008/05/01 

2008/08/13-2008/08/17 
2008/10/15-2008/10/19 
April- July 2008 



CAHA-3.5m/TWIN 
ING-WHT/ISIS 

CAHA-3.5m/TWIN 
ESO-NTT/EFOSC2 
ESO-VLT/FORS1 



Service 
P. M., S. 
S.B. 
H. H. 
A. T. 
Service 



G.. 



coverage is reached in the Northern hemisphere close to the 
galactic poles. SDSS data is widely used and therefore also well 
evaluated in terms of errors and accuracy (York et al. 2000: 
Abazajian et al. 2009 ). The SDSS data are supplemented by ad- 
ditional spectroscopic observations of appropriate quality from 
other sources. 



3.1. Colour selection and visual classification 

Hot subdwarfs are found most easily by applying a colour cut 
to Sloan photometry. All spectra of point sources with colours 
u - g < 0.4 and g — r < 0.1 were selected. This colour crite- 
rion corresponds to a limit in the Johnson photometric system of 
U - B < -0.57 (Jester et al. 120051 ). similar to the cut-off cho- 
sen by UV excess surveys, such as the Palomar Green survey 
(Green et al.|T986). The corresponding effective temperature of 



3 



S. Geier et al.: The MUCHFUSS project - Searching for hot subdwarf binaries with massive unseen companions 




g [mag] 



Fig. 3. Flux calibrated SDSS spectra of a single-lined sdB, a he- 
lium rich sdO and an sdB with main sequence companion visi- Fig 4 Heliocentric radial velocities of 1002 subdwarfs plotted 
ble in the spectrum. Note the different slopes of the sdB and the against g . magn i tude . The two dashed lines mark the RV cut of 
sdB+MS spectra. +100kms~ 1 



a BHB star is 15000 K (Castelli & Kurucz 2003), well be- 
low the observed range for sdB stars (> 20000 K). The limit of 
g - r - +0.1 corresponds to B - V = +0.3 (Jester et al. 120051 1. 
This ensures that sdBs in spectroscopic binaries are included if 
the dwarf companion is of spectral F or later, e.g. the sdB+F 
system PB 8783 at B - V = +0.13 and U - B = -0.65 (Koen 
et al. I1997K On the other hand the colour criteria exclude the 
huge number of QSOs (quasi stellar objects) which were the pri- 
ority objects of SDSS in the first place. We selected 48 267 point 
sources with spectra in this way. 

The spectra from SDSS are flux calibrated and cover the 
wavelength range from 3800 A to 9200 A with a resolution of 
R - 1800. Rebassa-Mansergas et al. (2007) verified the wave- 
length stability to be < 14.5 km s -1 from repeat sub-spectra us- 
ing SDSS observations of F-stars. We obtained the spectra of our 
targets from the SDSS Data Archive Servei0 and converted the 
wavelength scale from vacuum to air. The spectra were classified 
by visual inspection. 

In a first step extragalactic objects, spectra with low quality 
(S IN < 5) and unknown features have been excluded by visual 
inspection. In total we selected 10 811 spectra of 10 153 stars in 
this way. Fig. [2] (left panel) shows a two-colour plot of all se- 
lected objects. Classification was done by visual inspection of 
the spectra against reference spectra of hot subdwarfs and white 
dwarfs. Existence, width, and depth of helium and hydrogen ab- 
sorption lines as well as the flux distribution between 4000 and 
6000 A were used as criteria. Subdwarf B stars show broadened 
hydrogen Balmer and Hei lines, sdOB stars Hen lines in ad- 
dition, while the spectra of sdO stars are dominated by weak 
Balmer and strong He n lines depending on the He abundance. 
A flux excess in the red compared to the reference spectrum as 
well as the presence of spectral features such as the Mg i triplet 
at 5 170 A or the Ca n triplet at 8650 A were taken as indications 
of a late type companion (for a few examples see Fig. [3] for spec- 

3 das.sdss.org 



tral classification of hot subdwarf stars see the review by Heber 
[20091 . 

The sample contains 1100 hot subdwarfs in total. 725 be- 
long to the class of single-lined sdBs and sdOBs. Because dis- 
tinguising between these two subtypes from the spectral appear- 
ance alone can be difficult, we decided to treat them as one class. 
In addition we found 89 sdBs with cool companions. 198 stars 
are identified as single-lined sdOs, most of them enriched in he- 
lium. 9 sdOs have a main sequence companion. In 79 cases a 
unique classification was not possible. Most of these stars are 
considered as candidate sdBs of low temperature, which cannot 
be clearly distinguished from blue horizontal branch (BHB) stars 
or low mass white dwarfs of DA or DB type. 

Eisenstein et al. (2006) used a semi-automatic method for 
the spectral classification of white dwarfs and hot subdwarfs 
from the SDSS DR4. It is instructive to compare their sample 
to ours. Our colour cut-off is more restrictive and the confusion 
limit (S/N > 5) is brighter than that of Eisenstein et al. (2006i. 
Due to the redder colour cuts, blue horizontal branch stars en- 
ter the Eisenstein et al. sample, which we do not consider as hot 
subdwarf stars (see Heber 2009i. Applying our colour cuts to 
the hot subdwarf sample of Eisenstein et al. (12006) yields 691 
objects. The stars missing in our sample are mostly fainter than 
g = 19 mag as expected. Most recently, Kleinman (2010]) ex- 
tended the classifications to the SDSS DR7 and found 1409 hot 
subdwarf stars. Since no details are published, the sample can 
not be compared to ours yet. Considering our more restrictive 
colour cuts and confusion limit, the numbers compare very well 
with ours. This gives us confidence that our selection method is 
efficient. 

In Fig. [2] (right panel) only the subdwarf stars brighter than 
g = 18 mag are plotted. With less pollution by poor spectra, 
two sequences become clearly visible. The solid symbols mark 
single-lined sdBs and sdOs, while the open squares mark bi- 
naries with late type companions of most likely K and G type 
visible in the spectra. The contribution of the cool companions 
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Fig. 5. Radial velocity distribution of the hot subdwarf stars (see 
Fig- El- The bright sample (g < 16.5 mag, black histogram) con- 
tains a mixture of stars from the disk and the halo population. 
The faint sample {g > 16.5 mag, grey histogram) contains the 
halo population. The peak in the bright subsample around zero 
RV is caused by the thin disk population. The asymmetry in the 
faint subsample where negative RVs are more numerous than 
positive ones may be due to the presence of large structures in 
the halo and the movement of the solar system relative to the 
halo. 



shifts the colour of the star to the red. As can be seen in Figs. [2] 
the upper sequence also contains apparently single stars. Since 
the spectra are not corrected for interstellar reddening, some of 
these objects may show an excess in the red because of redden- 
ing rather than due to a cool companion. Amongst the faintest 
targets with noisy data, spectral features indicative of a late-type 
companion may have been missed as well as small excesses in 
the red. 

In Fig. [2] (right panel) we also compare the sample to syn- 
thetic colours suitable for hot subdwarf stars. We chose the grid 
of Castelli & Kurucz (20030 and selected models with high 
gravity (logg = 5.0). The models reproduce the lower envelope 
of the targets in the colour-colour-diagram very well for effective 
temperatures ranging from 20 000 to 50 000 K as expected for 
hot subdwarf stars. Different surface gravities, chemical com- 
positions and interstellar reddening are not accounted for, but 
would explain the observed scatter of the stars. 

It is interesting to note that there is an obvious lack of blue 
horizontal branch (BHB) stars with effective temperatures be- 
low 20 000 K compared to the sdBs with higher temperatures. 
This gap is hardly caused by selection effects, because the BHB 
stars are brighter than the sdBs in the optical. We conclude that 
the number density of BHB stars in the analysed temperature 
range must be much smaller than the one of sdBs. Newell (11973) 
was the first to report the existence of such a gap in the two- 
colour diagram of field blue halo stars, which was subsequently 
found to be also present in some globular clusters (Momany et 
al. 120041) . The reason for this gap remains unclear (see the review 
by Catelan 120091 . 

3.2. High radial velocity sample (HRV) 

The radial velocities of all identified hot subdwarf stars both 
single- and double-lined were measured by fitting a set of math- 
ematical functions (Gaussians, Lorentzians and polynomials) to 




4320 



4360 



http://wwwuser.oat.ts.astro.it/castelli/colors/sloan.html 



Fig. 6. H^-line of two consecutively taken individual SDSS spec- 
tra (Af = 0.056 d) of the sdB binary Jl 13840.68-003531.7. The 
shift in RV (^ 140km s ) between the two exposures is clearly 
visible. 



the hydrogen Balmer lines as well as helium lines if present 
using the FITSB2 routine (Napiwotzki et al. I2004ab and the 
Spectrum Plotting and Analysis Suite (SPAS) developed by H. 
Hirsch. Fig.|4]shows the RVs of 1002 hot subdwarf stars. 

Most of the known sdB binaries are bright objects (V 
10 - 14 mag) and the vast majority of them belongs to the 
Galactic disk population (Altmann et al. 2004). Due to the 
fact that these binary systems are close to the Sun they rotate 
around the Galactic centre with approximately the same veloc- 
ity. That is why the system velocities of most sdB binaries rela- 
tive to the Sun are low. One quarter of the known systems have 
\y\ < lOkms- 1 , 85% have \y\ < SOlons" 1 (see Table |A~TI 
Furthermore the RV semiamplitudes of these binaries are in most 
cases lower than lOOkms 1 (see Fig. [TJ. In order to filter out 
such normal thin-disk binaries we excluded sdBs with RVs lower 
than ±100kms _1 . 

Typical hot subdwarf stars fainter than g 17 mag have 
distances exceeding 4kpc and therefore likely belong to the 
Galactic halo population. Most of the stars in our sample are 
fainter than that (see Fig. @). The velocity distribution in the 
halo is roughly consistent with a Gaussian of 120 km s -1 disper- 
sion (Brown et al. 2005). Fig.[5]shows the velocity distribution 
of our sample dependent on the brightness of the objects. The 
distribution of the bright subsample (g < 16.5 mag) is roughly 
similar to the one of the faint subsample (g > 16.5 mag), the 
later extending to more extreme velocities and being somewhat 
asymmetric. Selecting objects with heliocentric radial velocities 
exceeding +100 km s _1 we aim at finding halo stars with extreme 
kinematics as well as close binaries with high RV amplitudes. 

Another selection criterion is the brightness of the stars. 
The accuracy of the RV measurements depends on the S/N of 
the spectra and the existence and strength of the spectral lines. 
Furthermore, the classification becomes more and more uncer- 
tain as soon as the S/N drops below ^ 10 and the probability 
of including DAs rises. Objects of uncertain type and RV (er- 
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Fig. 7. Probability for an sdB binary to host a massive compact 
companion and to be seen at sufficiently high inclination to un- 
ambiguously identify it from its binary mass function plotted 
against the RV shift within random times (solid curves, HRV 
sample) or on short timescales (dotted curve, RRV sample). 



rors larger than 50km s -1 ) have therefore been excluded. Most 
of the excluded objects are fainter than g = 19 mag. Altogether 
the target sample consists of 258 stars. 

Second epoch medium resolution spectroscopy was obtained 
starting in 2008 using ESO- VLT/FORS 1 (R 1800, /I = 
3730 - 5200 A), WHT/ISIS (R 4000,,} = 3440 - 5270 A), 
CAHA-3.5m/TWIN (R ^ 4000, A = 3460 - 5630 A) and ESO- 
NTT/EFOSC2 (R 2200, A = 4450 - 5 110 A). The journal of 
observations is given in Table Q] Up to now we have reobserved 
88 stars. We discovered 30 halo star candidates with constant 
high radial velocity (see Tillich et al. !2011b as well as 46 systems 
with radial velocities that were most likely variable. 

3.3. Rapid radial velocity variable sample (RRV) 

All SDSS spectra are co-added from at least three individual 
"sub-spectra" with typical exposure times of 15min. In most 
cases, the sub-spectra are taken consecutively; however, occa- 
sionally they may be split over several nights. In addition, some 
SDSS objects are observed more than once, either because the 
entire spectroscopic plate is re-observed, or because they are 
in the overlap area between adjacent spectroscopic plates. As 
a result, up to 30 sub-spectra are available for some objects. 
Hence, SDSS spectroscopy can be used to probe for radial veloc- 
ity variations, a method pioneered by Rebassa-Mansergas et al. 
(2007) to identify close white dwarf plus main-sequence bina- 
ries. We have obtained the sub-spectra for all sdBs brighter than 
g = 18.5 mag from the SDSS Data Archive Server. The quality 
of the individual spectra is not sufficient for our analysis in the 
case of even fainter stars. The object spectra were extracted from 
the FITS files for the blue and red spectrographs, and merged 
into a single spectrum using MIDAS. From the inspection of these 
data, we discovered 8 1 new candidate sdB binaries with radial 




RV [km s" 1 ] 

Fig. 8. Same as Fig. [7] except that the probability is plotted 
against RV at random time. 



velocity variations on short time scales, =2 0.02 - 0.07 d (see 
Fig. |6]for an example). 

The individual SDSS spectra are perfectly suited to search 
for close double degenerate binaries. Ongoing projects like 
SWARMS (Badenes et al. [20091 Mullally et al. 12009b fo cus on 
binaries with white dwarf primaries (see also Kilic et al. 2010; 
Marsh et al. 1201 Oi l and use a similar method. 

3.4. Selecting high mass companions 

Time resolved follow-up spectroscopy with a good phase cover- 
age is needed to determine the orbital solutions of the RV vari- 
able systems. In order to select the most promising targets for 
follow-up, we carried out numerical simulations and estimated 
the probability for a subdwarf binary with known RV shift to 
host a massive compact companion. We created a mock sample 
of sdBs with a close binary fraction of 50 %. 

We adopted the distribution of orbital periods of all known 
sdB binaries (see Table I A. Il l approximated by two Gaussians 
centered at 0.7 d (width 0.3 d) and 5.0 d (width 3.0 d) days 
and assumed that 82% of the binaries belong to the short pe- 
riod population. The short period Gaussian was truncated at 
0.05 d, which is considered the minimum period for an sdB bi- 
nary, because the subdwarf primary starts filling its Roche lobe 
for shorter periods and typical companion masses. Since sta- 
ble Roche lobe overflow and the accretion onto the compan- 
ion would dramatically change the spectra of these stars, we can 
safely presume that our sample does not contain such objects. 

The orbital inclination angles are assumed to be randomly 
distributed, but for geometrical reasons binaries at high incli- 
nations are more likely observed than binaries at low inclina- 
tions. To account for this, we used the method described in Gray 
( 1 19921 ) and adopted a realistic distribution of inclination angles. 

For the sdB mass the canonical value of 0.47 M Q was cho- 
sen. The distribution of companion masses was based on the re- 
sults by Geier et al. (2010b). The distribution of the low mass 
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companions was approximated by a Gaussian centered at 0.4 M 
(width 0.3 M ). The fraction of massive compact companions is 
estimated to 2 % of the close binary population based on binary 
population synthesis models (Geier et al.|2010b). The mass dis- 
tribution of these companions was approximated by a Gaussian 
centered at 2.0 M (width 1.0 M ). 

For the system velocities a Gaussian distribution with a dis- 
persion of 120 km s _1 typical for halo stars was adopted (Brown 
et al. 120051 ). Two RVs were taken from the model RV curves at 
random times and the RV difference was calculated for each of 
the 10 6 binaries in the simulation sample. This selection crite- 
rion corresponds to the HRV sample. For given RV difference 
and timespan between the measurements the fraction of systems 
with minimum companion masses exceeding 1 M was counted. 

In Fig. [7] the fraction of massive compact companions with 
unambiguous mass functions is plotted against the RV shift be- 
tween two measurements at random times (solid curve). It is 
quite obvious that binaries with high RV shifts are more likely to 
host massive companions. The probability for a high mass com- 
panion (> 1 M ) at high inclination is raised by a factor of ten as 
soon as the RV shift exceeds 200 km s . 

In order to check whether the selection of high velocities 
rather than high velocity shifts has an impact on the probability 
of finding sdB binaries with massive compact companions we 
used the same simulation. In Fig.[8]the fraction of these binaries 
is plotted against only one RV measurement taken at a random 
time. It can be clearly seen that the detection probability rises 
significantly for stars with high RVs. Selecting the fastest stars 
in the halo therefore makes sense when searching for massive 
compact companions to sdB. 

Since the individual SDSS spectra were taken within short 
timespans, another simulation was performed corresponding to 
the RRV sample. The first RV was taken at a random time, but 
the second one just 0.03 d later. The dotted curve in Fig.|7]illus- 
trates the outcome of this simulation. As soon as the RV shift 
exceeds 30 km s _1 within 0.03 d the probability that the compan- 
ion is massive rises to 10%. The reason why the probability 
does not increase significantly with increasing RV shift is that 
the most massive companions in our simulation have maximum 
RV shifts as high as 1000 km s~'. At the most likely periods of 
0.5 d the maximum RV shift within 0.03 d is then of the order 
of 100 km s -1 . Even higher RV shifts within short time are not 
physically plausible. 

Our simulation gives a quantitative estimate based on our 
current knowledge of the sdB binary populations. It has to be 
pointed out that these numbers should be considered as rough 
estimates at most. The observed period and companion mass 
distributions are especially affected by selection effects. The de- 
rived numbers are therefore only used to create a priority list and 
select the best targets for follow-up. 

3.5. Final target sample 

Our sample of promising targets consists of 69 objects in total. 
52 stars show significant RV shifts (> 30km s _1 ) within 0.02 - 
0.07 d and are selected from the RRV sample, while 17 stars 
show high RV shifts (100 - 300 kms" 1 ) within more than one 
day and are selected from the HRV sample (see Fig. [9j». 

In Geier et al. (1201 11 1 we showed that the SDSS spectra are 
well suited to determine atmospheric parameters by fitting syn- 
thetic line profiles to the hydrogen Balmer lines (Hp to H9) as 
well as He 1 and Hen lines. In order to maximize the quality of 
the data the single spectra were shifted to rest wavelength and 
coadded. The quality of the averaged spectra is quite inhomoge- 
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Fig. 9. Highest radial velocity shift between individual spectra 
plotted against time difference between the corresponding ob- 
serving epochs. The dashed horizontal line marks the selection 
criterion ARV > 100 kms 1 , the dotted vertical line the selection 
criterion Ar < 0. 1 d. All objects fulfilling at least one of these 
criteria lie outside the shaded area and belong to the top candi- 
date list for the follow-up campaign. The filled diamonds mark 
sdBs, while the blank squares mark He-sdOs. 

neous (S/N - 20 - 180, see Tabled, which affects the accuracy 
of the parameter determination. 

A quantitative spectral analysis was performed in the way 
described in Lisker et al. (2005 ) and Stroer et al. (2007 1. Due to 
the fact that our sample consists of different subdwarf classes, 
we used appropriate model grids in each case. For the hydrogen- 
rich and helium-poor (logy < -1.0) sdBs with effective tem- 
peratures below 30 000 K a grid of metal line blanketed LTE at- 
mospheres with solar metallicity was used. Helium-poor sdBs 
and sdOBs with temperatures ranging from 30 000 K to 40 000 K 
have been analysed using LTE models with enhanced metal line 
blanketing (O'Toole & Heber 2006). In the case of hydrogen- 
rich sdOBs with temperatures below 40 000 K showing moder- 
ate He-enrichment (logy = -1.0. .0.0) and hydrogen-rich sdOs 
metal-free NLTE models were used (Stroer et al. 2007 i. The 
He-sdOs have been analysed with NLTE models taking into ac- 
count the line-blanketing caused by nitrogen and carbon (Hirsch 
& Heber l2009l ). 

Spectral lines of hydrogen and helium were fitted by means 
of x 2 minimization using SPAS. The statistical errors have 
been calculated with a bootstrapping algorithm. Minimum er- 
rors reflecting systematic shifts when using different model grids 
(AJeff = 500 K; A logg = 0.05; A logy = 0.1, for a discus- 
sion see Geier et al. 2007 ) have been adopted in cases where the 
statistical errors were lower. Example fits for a typical sdB, an 
sdOB and a He-sdO star are shown in Fig.fTUl 

In addition to statistical uncertainities systematic effects have 
to be taken into account in particular for sdB stars. The higher 
Balmer lines (H e and higher) at the blue end of the spectral 
range are very sensitive to changes in the atmospheric param- 
eters. However, the SDSS spectral range restricts our analysis to 
the Balmer lines from Hp to H9. In high S/N data these lines are 
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Fig. 10. Example fits of hydrogen and helium lines with model spectra for an sdB (left panel), an sdOB (middle panel) and a He-sdO 
star (right panel). The atmospheric parameters of these stars are given in Tables[3]and|4] 



sufficient to measure accurate parameters as has been shown in 
Geier et al. (1201 U . In spectra of lower quality the bluest lines 
(H9 and Hs) are dominated by noise and cannot be used any 
more. In order to check whether this leads to systematic shifts in 
the parameters as reported in Geier et al. (12010b) we made use 
of the individual SDSS spectra. We chose objects with multiple 
spectra, which have an S/N comparable to the lowest quality data 
in our sample (^ 20). The atmospheric parameters were obtained 
from each individual spectrum. Average values of T e s and log g 
were calculated and compared to the atmospheric parameters de- 
rived from the analysis of the appropriate coadded spectrum. For 
effective temperatures ranging from 27 000 K and 39 000 K no 
significant systematic shifts were found. This means that the er- 
ror is dominated by statistical noise. However, for temperatures 
as low as 25 000 K systematic shifts of the order of -2500 K 
in r e ff and -0.35 in logg are present. For sdBs with low effec- 
tive temperatures and signal-to-noise, the atmospheric parame- 
ters are therefore systematically underestimated. Only three stars 
in our sample have temperatures in this range. Since their coad- 
ded spectra are of reasonable quality (S/N = 34 - 167), sys- 
tematic shifts should be negligible in these cases. Because all 
important lines of Hei and Hen are well covered by the SDSS 
spectral range, systematic effects should be negligible in the case 
of He-rich sdO/Bs as well. 

The parameters of the sample are given in Tables|3] and 
Seven stars have already been analysed in Geier et al. (1201 1) . 
The sample consists of 38 hydrogen rich sdBs, 13 sdOBs and 
3 hydrogen rich sdOs. Thirteen stars are helium rich sdOs (He- 
sdOs) and J134352.14+394008.3 belongs to the rare class of he- 
lium rich sdBs. 

Our SDSS sample reaches down to fainter magnitudes and 
hence, larger distances than any previous survey. In an ongoing 
project Green et al. (20081 analyse all hot subdwarfs from the 
PG survey down to ^ 14.0 mag. The sample of hot subdwarf 
stars analysed in the course of the SPY survey reaches down to 
^ 16.5 mag (Lisker et al. 2005; Stroer et al. 2007 1, quite similar 
to the sample of sdBs from the Hamburg Quasar Survey analysed 
by Edelmann et al. (120031 . 

Spectroscopic distances to our stars have been calculated 
as described in Ramspeck et al. (2001) assuming the canonical 
mass of 0.47 M for the subdwarfs and using the formula given 



by LuptorQ to convert SDSS-g and r magnitudes to Johnson 
V magnitudes. Again interstellar reddening has been neglected. 
The distances range from 1 kpc to > 16kpc. Since the SDSS 
footprint is roughly perpendicular to the Galactic disk, these dis- 
tances tell us something about the population membership of our 
stars. These subdwarfs most likely belong to the thick disk or the 
halo with small contributions of thin disk stars. 

Fig.fTTIshows a T e s - \ogg diagram of the top target sample. 
Most of our stars were born in an environment of low metallicity 
(thick disk or halo). Dorman et al. ( 1993 1 calculated evolutionary 
tracks for different metallicities of the subdwarf progenitor stars. 
For lower metallicities, the evolutionary tracks and with them the 
location of the EHB, are shifted towards higher temperatures and 
lower surface gravities. In Fig.QT]the T e ff - logg diagram is su- 
perimposed with evolutionary tracks and an EHB calculated for 
a subsolar metallicity of log z = -1 .48, which is consistent with 
a mixture between thick disk and halo population. Evolutionary 
tracks for solar metallicity are given in Fig. [12]for comparison. 

Most of the sdB stars with hydrogen-rich atmospheres are 
found on or slightly above the EHB band implying an evolution- 
ary status as core helium-burning EHB or shell helium-burning 
post-EHB stars. The sample contains only three hydrogen rich 
sdOs, which are thought to be evolved post-EHB stars in a tran- 
sition state. The He-sdOs cluster near the HeMS at temperatures 
of 45 000 K. This is fully consistent with the results from the 
PG and the SPY s urveys (Green et al. 120081 Lisker et al. 120051 
Stroer et al. 12007) and illustrates that our sample is not biased 
(see Fig. IT2T>. 

Compared to other studies we find only a few stars with tem- 
peratures lower than 27 000 K. Furthermore, the scatter around 
the EHB seems to be systematically shifted towards higher tem- 
peratures and lower surface gravities. According to our study of 
systematic errors in the parameter determination, it is unlikely 
that this causes the effect. However, higher quality data would 
be necessary to verify this. Another possible explanation might 
be related to the volume of the sample. Since hot subdwarfs of 
lower temperature are brighter in the optical range because of 
the lower bolometric correction, we may already see all of them 
in a fixed volume, while the fraction of hot stars is still rising at 
fainter magnitudes. 



5 http://www.sdss.org/dr6/algorithms/sdssUBVRITransform.html 
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Fig. 11. r e jf - logg diagram of our target sample. The helium 
main sequence (HeMS) and the EHB band (limited by the zero- 
age EHB, ZAEHB, and the terminal-age EHB, TAEHB) are su- 
perimposed with EHB evolutionary tracks for subsolar metallic - 
ity (logz = -1.48) from Dorman et al. ([1993). 

In Fig. Q~3] the helium abundance is plotted against effective 
temperature. The general correlation of helium abundance with 
effective temperature and the large scatter in the region of the 
sdB stars have been observed in previous studies as well. Two 
sequences of helium abundance among the sdB stars as reported 
by Edelmann et al. (2003 ) could not be identified. 

One has to keep in mind that our sample consists of RV vari- 
able stars only. In Fig. QT|a lack of such stars at the hot end of 
the EHB is visible. Green et al. (2008) reported similar system- 
atics in their bright PG sample. The reason for this behaviour 
is not fully understood yet. According to the model of Han et 
al. (120021 120031 and Han d2008b sdBs with thin hydrogen en- 
velopes situated at the hot end of the EHB may be formed after 
the merger of two helium WDs. Since merger remnants are sin- 
gle stars, they are not RV variable. 

The top target sample includes 13 He-sdOs where RV shifts 
of up to lOOkms -1 have been detected within short timespans 
of 0.01 - 0.1 d. In total 20 He-sdOs show signs of RV variabil- 
ity. This fraction was unexpected since the fraction of close bi- 
nary He-sdOs from the SPY sample turned out to be 4% at most 
(Napiwotzki|2008ll|!| 

4. Summary and Outlook 

In this paper we introduced the MUCHFUSS project, which 
aims at finding sdBs in close binaries with massive compact 
companions. We identified 1100 hot subdwarf stars from the 
SDSS by colour selection and visual inspection of their spectra. 
Stars with high absolute radial velocities have been selected to 
efficiently remove normal sdB binaries from the thin disk popu- 
lation and have been reobserved. 46 binary candidates with sig- 

6 Green et al. 02008 ) suggested that the binary fraction of He-sdO 
stars may be comparable to the binary fraction of sdBs. 



Fig. 12. T e ff - log g diagram of the hot subdwarfs from the SPY 
project (Lisker et al. 2005, Stroer et al. 2007i. The helium main 
sequence (HeMS) and the EHB band (limited by the zero-age 
EHB, ZAEHB, and the terminal-age EHB, TAEHB) are super- 
imposed with EHB evolutionary tracks for solar metallicity from 
Dorman et al. JT9931 . 

nificant RV shifts have been found. From the analysis of indi- 
vidual SDSS spectra, 81 additional stars with RV shifts on short 
timescales have been found. 

Targets for follow-up spectroscopy have been selected using 
numerical simulations based on the properties of the known sdB 
close binary population and theoretical predictions about the rel- 
ative fraction of massive compact companions. 69 binaries with 
high RV shifts as well as significant RV shifts on short timescales 
have been selected as good candidates for massive compact com- 
panions. Atmospheric parameters, spectroscopic distances and 
population memberships have been determined. 

The multi-site follow-up campaign started in 2009 and is be- 
ing conducted with medium resolution spectrographs mounted 
at several different telescopes of mostly 4-m-class. First results 
are presented in Geier et al. (1201 1) . 
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Fig. 13. Helium abundance logy plotted against effective tem- 
perature (see Tables 1314) . The solid horizontal line marks the 
solar value. Lower and upper limits are marked with upward and 
downward triangles. 
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Table 2. Priority targets for follow-up. Besides the names, the g magnitudes, the number of individual spectra and the S/N of the 
coadded spectra at ^ 4100 A are given. 



Object 




g 


No. 


S/N 


Object 




g 


No. 


S/N 


J002323.99-002953.2 


PB5916 


15.3 


16 


116 


J150513.52+110836.6 


PG 1502+113 


15.1 


4 


90 


J012022.94+395059.4 


FBS 01 17+396 


15.2 


8 


100 


J 150829.02+494050.9 




17.3 


3 


50 


J012739.35+404357.8 




16.5 


8 


59 


J151415.66-012925.2 




16.8 


5 


48 


J052544.93+630726.0 




17.6 


3 


35 


J152222.15-013018.3 




17.7 


5 


28 


J074534.16+372718.5 




17.6 


5 


26 


J152705.03+1 10843.9 




17.1 


5 


39 


J075937. 15+541022.2 




17.5 


3 


27 


J153411. 10+543345.2 


WD 1532+547 


16.7 


8 


52 


J082053.53+000843.4 




14.9 


6 


103 


J155628.34+01 1335.0 




16.0 


8 


92 


J083006. 17+475 150.4 




15.8 


5 


95 


J161 140.50+201857.0 




18.2 


5 


20 


J085727.65+424215.4 


US 1993 


18.3 


4 


21 


J161817.65+120159.6 




17.8 


4 


18 


J092520.70+470330.6 




17.4 


3 


33 


J 162256. 66+47305 1.1 




16.0 


4 


72 


J094856.95+334151.0 


KUV 09460+3356 


17.4 


3 


46 


J163702.78-01 1351.7 




17.1 


12 


46 


J095229.62+301553.6 




18.2 


3 


20 


J164326.04+330113.1 


PG 1641+331 


16.1 


3 


55 


J095238.93+252159.7 




14.5 


4 


113 


J165404.26+303701.8 


PG 1652+307 


15.1 


4 


167 


J100535.76+223952.1 




18.1 


4 


28 


J170645.57+243208.6 




17.5 


3 


39 


J102151.64+301011.9 




18.0 


12 


34 


J170810.97+244341.6 




18.2 


3 


16 


J103549.68+092551.9 




16.0 


3 


59 


J171617.33+553446.7 


SBSS 1715+556 


16.9 


8 


39 


Jl 10215.97+521858.1 




17.2 


3 


44 


J171629.92+575121.2 




17.9 


4 


21 


Jl 10445.01+092530.9 




16.0 


4 


40 


J172624. 10+274419.3 


PG 1724+278 


15.7 


4 


107 


Jl 12242.69+613758.5 


PG 11 19+619 


15.1 


3 


87 


J1745 16.32+244348.3 




17.4 


3 


22 


Jl 12414.45+402637.1 




17.7 


3 


21 


J175 125.67+255003.5 




17.2 


4 


50 


Jl 13303.70+290223.0 




17.4 


3 


34 


J202313.83+131254.9 




17.0 


3 


33 


Jl 13418.00+015322.1 


LBQS 1131+0209 


17.7 


6 


30 


J202758.63+773924.5 




17.7 


3 


22 


Jl 13840.68-003531.7 


PG 1136-003 


14.2 


10 


174 


J204300.90+002145.0 




17.6 


9 


50 


Jl 13935.45+614953.9 


FBS 1136+621 


16.8 


3 


34 


J204448.63+ 153638.8 




17.7 


7 


50 


J115358.81+353929.0 


FBS 1151+359 


16.3 


3 


48 


J204546.8 1-054355.6 




17.8 


4 


29 


Jl 15716.37+612410.7 


FBS 1154+617 


16.9 


5 


34 


J204613. 40-045418. 7 




16.0 


3 


120 


J125702.30+435245.8 




17.9 


3 


18 


J204940.85+ 165003.6 




17.7 


7 


35 


J 130059.20+0057 11.7 


PG 1258+012 


16.3 


3 


47 


J210454.89+110645.5 




17.2 


4 


37 


J130439.57+312904.8 


LB 28 


16.8 


3 


42 


J21 1651.96+003328.5 




17.7 


3 


19 


J133638.81+111949.4 




17.0 


3 


32 


J215648.71+003620.7 


PB 5010 


17.7 


3 


22 


J134352.14+394008.3 




18.1 


3 


19 


J225638.34+065651.1 


PG 2254+067 


15.1 


3 


86 


J135807.96+261215.5 




17.7 


4 


23 


J232757.46+483755.2 




15.6 


3 


92 


J140545.25+014419.0 


PG 1403+019 


15.6 


3 


81 


J233406.1 1+462249.3 




17.4 


3 


35 


J141549.05+111213.9 




15.8 


3 


82 


J234528.85+393505.2 




17.3 


3 


37 


J143 153.05-002824.3 


LBQS 1429-0015 


17.8 


8 


34 













Table 3. Priority targets for follow-up (HRV subsample). fThe binary system has been analysed in Geier et al. (1201 11 1. 



Object 






Class 






l°gg 






logy 


d 


ARV 




At 










[K] 














[kpc] 


[kms~ 


'] 


[d] 


J102151 


.64+301011 


.9 


sdB 


30700 


±500 


5.71 


± 


0. 


06 


< -3.0 




277 ± 


51 


14.936 


J150829 


.02+494050 


.9 


sdB 


28200 


±600 


5.34 


± 


0. 


09 


-2.0 ±0.2 


6.4+»; 


211 ± 


18 


2161.429 


J095229 


.62+301553 


.6 


sdB 


35200 


± 1200 


5.05 


± 


0. 


17 


< -3.0 


16.0™ 


198 ± 


40 


1155.766 


Jl 13840 


.68-003531 


•7t 


sdB 


30800 


±500 


5.50 


± 


0. 


09 


-3.0 ±0.2 


1 3 +0 - 2 ' 

-0.1 


182 ± 


12 


0.973 


J 165404 


.26+303701 


.8t 


sdB 


24400 


±800 


5.32 


± 


0. 


.11 


-2.3 ±0.3 


1.9+03 
4 8 +1 ^ 


181 ± 


9 


1795.144 


J 152222 


.15-013018 


.3 


sdB 


24800 


± 1000 


5.52 


± 


0. 


15 


-2.6 ±0.5 


173 ± 


36 


3.001 


J150513 


.52+110836 


.6t 


sdB 


33300 


±500 


5.80 


± 


0. 


10 


-2.4 ±0.3 


15® 


154 ± 


12 


0.957 


J002323 


.99-002953 


.2t 


sdB 


30100 


±500 


5.62 


± 


0. 


08 


-2.2 ±0.2 


1 8 +0 - 2 


130 ± 


14 


82.784 


J202313 


.83+131254 


.9 


sdB 


29600 


±600 


5.64 


± 


0. 


14 


-2.1 ±0.1 




124 ± 


21 


1202.795 


JO 12022 


.94+395059 


.4 


sdB 


28900 


±500 


5.51 


± 


0. 


08 


-3.0 ±0.4 


3s 


114 ± 


11 


360.973 


J202758 


.63+773924 


.5 


sdO 


46200 


±3200 


5.48 


± 


0. 


18 


-2.8 ±0.9 




114 ± 


48 


1.960 


J095238 


.93+252159 


.7 


sdB 


27800 


±500 


5.61 


± 


0. 


08 


-2.64 ±0.1 




111 ± 


10 


2.918 


J161140 


.50+201857 


.0 


sdOB 


36900 


±700 


5.89 


± 


0. 


13 


-1.2 ±0.1 




108 ± 


36 


0.947 


J 164326 


.04+330113 


.1 


sdB 


27900 


±500 


5.62 


± 


0. 


07 


-2.3 ±0.2 


2 4 +0 - 2 


108 ± 


11 


1.990 


J204448 


.63+153638 


.8 


sdB 


29600 


±600 


5.57 


± 


0. 


09 


-2.2 ±0.1 


101 ± 


19 


3.049 


J083006 


.17+475150 


.4 


sdB 


25200 


±500 


5.30 


± 


0. 


05 


-3.3 ±0.7 


2-8-o:l 


95 ± 14 


3.961 


J204940 


.85+165003 


.6 


He-sdO 


43000 


±700 


5.71 


± 


0. 


13 


> +2.0 




85 ± 19 


5.932 
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Table 4. Priority targets for follow-up (RRV subsample). fThe binary system has been analysed in Geier et al. (1201 U . ^Atmospheric 
parameters (T e g = 39400 K, logg = 5.64, logy = -0.55) have been determined by Stroer et al. (12007V 



J085727 


,65+424215 


.4 


He-sdO 


39500 


± 


1900 


5 


,63 


±0. 


24 


+0.2 ± 0.2 


8. 


.7 - 9 


111 : 


+ 46 


0.066 


J161817 


.65+120159 


.6 


sdB 


32100 


± 


1000 


5 


.35 


±0. 


23 


— 


8. 


1+2.8 
' 1 -2.1 

O+U.J 
■ -0.3 


105: 


+ 31 


0.043 


J232757 


,46+483755 


.2 


He-sdO 


64700 


± 


2000 


5 


,40 


±0. 


08 


> +2.0 


4. 


105: 


+ 24 


0.016 


J 162256 


,66+473051 


.1 


sdB 


28600 


+ 


500 


5 


,70 


±0. 


11 


-1.81 +0.1 


2. 


101 : 


+ 15 


0.037 


J 163702 


,78-011351 


.7 


He-sdO 


46100 


± 


700 


5 


,92 


±0. 


22 


> +2.0 


3. 


Q + l . 1 

■ 8 -0.9 


101 : 


+ 55 


0.085 


Jl 13303 


,70+290223 


.0 


sdB/DA 


— 






— 








— 


— 


95 ± 


35 


0.016 


J 135 807 


,96+261215 


.5 


sdB 


33500 


± 


600 


5 


,66 


±0. 


10 


< -2.0 


5. 




87 ±29 


0.030 


Jl 12242 


,69+613758 


.5 


sdB 


29300 


+ 


500 


5 


.69 


± 0. 


10 


-2.3 + 0.3 


1 


83 ± 


20 


0.047 


J153411 


,10+543345 


.2 


sdOB 


34800 


± 


700 


5. 


.64 


± 0. 


09 


-2.6 ± 0.3 


3, 


> 


83 ± 


29 


0.018 


J082053 


,53+000843 


.4 


sdB 


26700 


± 


900 


5 


.48 


± 0. 


10 


-2.0 ± 0.09 


1. 


77 ± 


11 


0.047 


J170810. 


,97+244341 


.6 


sdOB 


35600 


± 


800 


5, 


.58 


±0. 


14 


-0.8 ± 0.1 


8, 


+i'S 


76 ± 


33 


0.013 


J094856 


,95+334151 


.0 


He-sdO 


51000 


± 


1200 


5, 


.87 


± 0. 


12 


+ 1.8 ± 0.5 


5. 




75 ± 


17 


0.012 


J204613 


,40-045418 


■7t 


sdB 


31600 


± 


600 


5 


.55 


± 0. 


10 


-3.7 ± 0.6 


2. 


hi 


70 ± 


13 


0.030 


J215648. 


,71+003620 


.7 


sdB 


30800 


± 


800 


5, 


.77 


± 0. 


12 


-2.2 ± 0.3 


4. 




69 ± 


21 


0.011 


J074534 


.16+372718 


.5 


sdB 


37500 


± 


500 


5 


.90 


± 0. 


09 


< -3.0 


4. 


6 


65 ± 


19 


0.036 


J143153. 


,05-002824 


.3 


sdOB 


37300 


± 


800 


6 


.02 


± 0. 


16 


-0.8 ± 0.1 


4. 


j+u.y 
-Q.8 


65 ± 


22 


0.012 


J171629. 


,92+575121 


.2 


sdOB 


35400 


± 


1000 


5, 


.60 


± 0. 


18 


-0.7 ± 0.1 


7. 


fjl 

Q+i .y 

-14 

Q+1.3 

-2.2 


65 ± 


16 


0.013 


Jl 12414 


.45+402637 


.1 


He-sdO 


47100 


± 


1000 


5 


.81 


± 0. 


23 


+ 1.7 ± 0.7 


5. 


63 ± 


22 


0.021 


J 125702. 


,30+435245 


.8 


sdB 


28000 


± 


1100 


5, 


.77 


± 0. 


17 


< -3.0 


4. 


63 ± 


28 


0.010 


Jl 10215 


,97+521858 


.1 


He-sdO 


56600 


± 


4200 


5, 


.36 


± 0. 


22 


> +2.0 


8. 


62 ± 


11 


0.033 


J151415 


,66-012925 


.2 


He-sdO 


48200 


± 


500 


5 


.85 


± 0. 


08 


+ 1.7 ± 0.4 


3. 


6 +0 - 4 

u -0.3 


62 ± 


22 


0.016 


J204300. 


,90+002145 


.0 


sdO 


40200 


± 


700 


6 


.15 


± 0. 


13 


-1.3 ± 0.4 


3. 


(.+0.6 

9^ 


61 ± 


13 


0.016 


J171617 


,33+553446 


.7 


sdB 


32900 


± 


900 


5 


.48 


± 0. 


09 


< -3.0 


4. 


60 ± 


24 


0.048 


J2 10454. 


,89+110645 


.5 


sdOB 


37800 


± 


700 


5, 


.63 


± 0. 


10 


-2.4 ± 0.2 


4. 




58 ± 


19 


0.023 


Jl 15358 


,81+353929 


.0 


sdOB 


29400 


± 


500 


5 


.49 


± 0. 


06 


-2.5 ± 0.3 


3. 




56 ± 


12 


0.022 


J174516 


,32+244348 


.3 


He-sdO 


43400 


± 


1000 


5 


.62 


± 0. 


21 


> +2.0 


6. 


+ 18 
-1.4 


55 ± 


28 


0.016 


J134352. 


,14+394008 


.3 


He-sdB 


36000 


± 


2100 


4 


.78 


± 0. 


30 


-0.2 ± 0.2 


8. 




52 ± 


34 


0.022 


Jl 15716 


,37+612410 


.7 


sdB 


29900 


± 


500 


5 


.59 


± 0. 


08 


-3.2 ± 0.8 


4. 


r\+0.5 
"-0.4 


51 ± 


34 


0.049 


J133638. 


,81 + 111949 


.4 


sdB 


27500 


± 


500 


5, 


.49 


± 0. 


08 


-2.7 +0.2 


4. 




48 ± 


17 


0.030 


J211651 


,96+003328 


.5 


sdB 


27900 


± 


800 


5 


.78 


± 0. 


15 


-3.9 ± 0.7 


4. 


^+0'6 


48 ± 


23 


0.016 


J 170645. 


,57+243208 


.6 


sdB 


32000 


± 


500 


5, 


.59 


± 0. 


07 


< -4.0 


5. 


'^-0.5 


46 ± 


14 


0.013 


J175125. 


,67+255003 


.5 


sdB 


30600 


± 


500 


5 


.48 


± 0. 


08 


< -3.8 


5. 


r\+0.6 

tlii 


46 ± 


14 


0.034 


J012739. 


,35+404357 


.8 


sdO 


48300 


± 


3200 


5, 


.67 


± 0. 


10 


-1.3 +0.2 


4. 


45 ± 


17 


0.037 


Jl 13418 


,00+015322 


.1 


sdB 


29700 


± 


1200 


4 


.83 


± 0. 


16 


< -4.0 


1. 




45 ± 


24 


0.076 


J172624. 


,10+274419 


,3f 


sdOB 


33500 


± 


500 


5, 


.71 


± 0. 


09 


-2.2 ± 0.1 


2. 


"^-0.2 


45 ± 


16 


0.047 


J155628. 


,34+011335 


.0 


sdB 


32700 


± 


600 


5, 


.51 


± 0. 


08 


-2.9 + 0.2 


3. 




44 ± 


15 


0.068 


J103549 


,68+092551 


.9 


He-sdO 


48100 


± 


600 


6 


.02 


± 0. 


13 


> +2.0 


2. 




43 ± 


12 


0.021 


J141549. 


,05+111213 


.9 


He-sdO 


43100 


± 


800 


5 


.81 


± 0. 


17 


> +2.0 


2. 


4 ;8:1 


43 ±7 


0.023 


J152705 


,03+110843 


.9 


sdOB 


37600 


± 


500 


5 


.62 


± 0. 


10 


-0.5 ± 0.1 


4. 


,8 ;8-6 


43 ± 


14 


0.054 


J052544 


,93+630726 


.0 


sdOB 


35600 


± 


800 


5, 


.85 


± 0. 


10 


-1.6 ± 0.2 


4. 




42 ± 


17 


0.026 


J100535 


.76+223952 


.1 


sdB 


29000 


± 


700 


5 


.43 


± 0. 


13 


-2.7 ± 0.2 


7. 




41 ± 


18 


0.019 


J204546. 


,81-054355 


.6 


sdB 


35500 


± 


500 


5, 


.47 


± 0. 


09 


-1.4 ±0.2 


7. 




41 ± 


18 


0.013 


J092520. 


,70+470330 


.6 


sdB 


28100 


± 


900 


5, 


.17 


± 0. 


15 


-2.5+0.2 


7. 


-1.4 


40 ± 


13 


0.012 


J075937 


,15+541022 


.2 


sdB 


31300 


± 


700 


5 


.30 


±0. 


10 


-3.3 ±0.3 


7, 


5+t.i 


38 ± 


13 


0.012 


J234528. 


,85+393505 


.2 


He-sdO 


47900 


± 


800 


6 


.07 


±0. 


14 


> +2.0 


3, 


c+0.6 
J -0.5 


37 ± 


14 


0.012 


J130439 


,57+312904 


.8 


sdOB 


38100 


± 


600 


5 


.69 


±0. 


12 


-0.4 ±0.1 


4. 


i +0.6 


36 ± 


12 


0.037 


J 130059 


,20+005711 


.7* 


He-sdO 


40700 


± 


500 


5, 


.53 


±0. 


10 


-0.6 ±0.1 


3, 


9+0.5 

o+8i 
' S -0.3 


36 ± 


16 


0.012 


Jl 10445 


,01+092530 


.9 


sdOB 


35900 


± 


800 


5 


.41 


±0. 


07 


-2.1 ±0.4 


3. 


34 ± 


14 


0.040 


Jl 13935 


,45+614953 


.9 


sdB 


28800 


± 


900 


5, 


.27 


±0. 


15 


-2.8 ± 0.3 


4. 


q+Il 
■ y -a.9 


31 ± 14 


0.011 


J233406 


,11+462249 


.3 


sdOB 


34600 


± 


500 


5 


.71 


±0. 


09 


-1.3 ±0.1 


4. 


q+o.6 


31 ± 


14 


0.025 


J225638 


,34+065651 


■ It 


sdB 


28900 


± 


600 


5 


.58 


±0. 


11 


-3.0 ±0.2 


1. 


6^ 


27 ± 


11 


0.031 


J 140545. 


,25+014419 


.0 


sdB 


27300 


± 


800 


5, 


.37 


±0. 


16 


-1.9 ±0.2 


2. 


c+0.6 
J -0.5 


25 ± 


10 


0.026 
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Appendix A: Close binary subdwarfs from literature 
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S. Geier et al.: The MUCHFUSS project - Searching for hot subdwarf binaries with massive unseen companions 

Table A.l. Orbital parameters of all known hot subdwarf binaries from literature. The superscript p denotes sdB pulsators, r binaries 
where with reflection effect, ec eclipsing systems and el systems with light variations caused by ellipsoidal deformation. fPost-RGB 
stars without core helium-burning. ^Double-lined binary consisting of two helium rich sdBs. The RV semi-amplitudes of both 
components are given. 



\J UJ CC I 


p 

r 


7 




Reference 




Ml 


[Km s J 


[Km s J 




PG 0850+ 170 


27.815 


32.2 ±2.8 


33.5 ±3.3 


Morales-Rueda et al. 2003a 


PG 1619+522 


15.3578 


-52.5 ± 1.1 


35.2 ± 1.1 


Morales-Rueda et al. 2003a 


PG 11 10+294 


9.4152 


-15.2 + 0.9 


58.7 ± 1.2 


Morales-Rueda et al.l2003al 


Feige 108 


8.7465 


45.8 ± 0.6 


50.2 ± 1.0 


Edelmann et al. 120041 


PG 0940+068 


8.330 


-16.7 ± 1.4 


61.2 ± 1.4 


Maxtedetal.l2000bl 


PHL861 


7.44 


-26.5 ± 0.4 


47.9 + 0.4 


Karl et al.l2006l 


HE 1448-0510 


7.159 


-45.5 + 0.8 


53.7 ± 1.1 


Karl etal. 120061 
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